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In this paper, a hybrid power system is proposed to eliminate any interruption 
caused to absence of power supply network by the Fusion between the wind 
source and the photovoltaic source to create a hybrid system with internal 
power supply autonomy. The hybrid power system includes a solar and wind 
source whose wind energy conversion system is based on a doubly fed 
induction generator. The photovoltaic conversion system extracts the 
maximum power from the sun. A part of this power is used to supply the 
rotor of a doubly fed induction generator; the rest is injected in the grid. The 
connection between the two conversion systems is made by means of static 
converters (chopper/inverter). In the photovoltaic conversion chain, a 
controller for tracking the maximum power point is designed using the direct 
search approach the Perturb and observe method. The wind energy 
conversion system contains a doubly fed asynchronous generator with a 
control system of generated powers by the technique of the fuzzy logic. The 


simulation was done to prove the validity of the strategy and of the control 
method used for the control of the power generated by the hybrid system. 
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1. INTRODUCTION 

Renewable energies have become very promising sources of future energy resources in the world, 
primarily because of their own nature, and also because of the diminishing reserves of traditional resources 
(oil, gas, etc.). .), Demographic and economic growth [1, 2]. Where the use of these sources to generate 
electricity has been seen since the last decade [3]. The main advantage of renewable sources is the absence of 
harmful emissions. One way to generate electricity from renewable sources is to use wind and photovoltaic 
power, mainly because it is clean [4]. The hybrid system is achieved in renewable energies consisting of 
wind energy and photovoltaic (PV) energy for parallel performance of these two sources, ensuring DFIG 
training and maximizing power generation with better efficiency [5, 6] through the integration of foggy logic 
technology to ensure control of powers Deposited towards the network. 

This control technology is simple in its implementation and its relatively satisfactory efficiency (for 
some applications), and the ease of operation and low cost of these controllers are the main reasons for its 
success [7]. Fuzzy logic is therefore one of the important branches of artificial intelligence among all 
intelligent techniques, the simplest method of control to be implemented in practice. [8, 9]. 
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The work presented in this paper relates to the implementation of the DFIG wind turbine control 
strategy and the supply of rotor from a PV generator. This system is controlled by foggy logic technology. 
The regulation of this paper is as follows: In Section 2, we introduce the proposed mixed system (PV/wind). 
3rdsection was designed to model hybrid system components. In Section 4, a control strategy is presented. 
In Section 5, the principle of precise logical control is described and used in the design of misty logical 
controllers. The photovoltaic system model is achieved in Section VI. In Section 7, we describe the MPPT 
structure applied to the PV system. The simulation results and conclusions are summarized in the last 
sections 


2. PRESENTATION OF THE HYBRID SYSTEM (PHOTOVOLTAIC/WIND) 

The hybrid power system contains two conversion chains. The first chain is a wind power 
conversion system, containing a DFIG connected directly to the networks. The second chain is a photovoltaic 
conversion system. The power generated by the DFIG is controlled by the technique of the fuzzy logic. 
The PV generator is used to supply the DFIG rotor, while the rest of the generated power is injected into the 
grid. Figure 1 shows the block diagram of the proposed hybrid power system. 


Electrical networks 








Figure 1. Block diagram the hybrid power system 


3. MODELING OF TNE CONVERSION HYBRID SYSTEM 
3.1. Model of the turbine 
The aeroturbine is composed of three pales of length R on one tree, a turbine and a gearbox with a 
gain G. The generator shaft is modeled by the following equation [10]: 





d 
Jo Qmec = C g~Cem~ Smee (1) 
P 
Cg= Caer C4, = G2 mec = GQur (2) 
tur 


Where: J is the total inertia of rotating parts (kg * m7); f is the viscous friction coefficient; C.,, is the 
electromagnetic torque of the generator (N * m); C, is the mechanical torque in the output of the multiplier; 
Caer is the aerodynamic torque of the turbine; G is the multiplier gain; Qnec is generator speed; Qy, is the 
speed of the turbine. In the wind turbine model, the aerodynamic power output is given as follows: 
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Where p is the air density (Kg/m3), v is the wind speed (m/s), and C, is the power coefficient. 
The power coefficient C, presents the portion of wind power captured by the turbine. It depends on 
two parameters, namely tip speed ratio 4 and pitch angle of the blades f: [11, 12]: 


m(A+0.1) 


C= 05 -o00167(a-2)in| FE) 


|-noniss(a—y(0-2) (4) 


Where: 4 is the speed ratio is defined as the ratio between the linear speed of the blades and the wind speed; 
f is the orientation angle of the blades. 
3.2. Modeling of DFIG 


The Park model of the machine is described by the following equations [13, 14]: The electric 
equations: 


. d 
Vds = Rsidst “eds - Os gs 


: d 
Vas — R sigs + “eds m Os as 





d (5) 
Vdr — Rriar+ 7 $ar-(@s- Or) b gy 
d d 
Var = Rrigr+ Fb grt (Os- Or) bay 
The flux equations: 
P ds = L slds a M srldr 
P as = L sigs +M srlgr is 
P ay =Lyidr+M srids 
P ar ~ L rigr +M srlqs 
The electromagnetic torque equation: 
M sr : ’ 
Cem = Pa (bgsiar Paster] (7) 
Ss 


where: R; is the stator resistance, R, is the rotor resistance, Ls is the stator inductance, Ms is mutual 
inductance, ¢as, @gs are respectively the direct and quadrature fluxes of the stator, ¢a,, dg are respectively the 
direct and quadrature fluxes of the rotor, ia, igs are respectively direct and quadrature stator currents, ia, igr 
are respectively direct and quadrature rotor currents, p is number of even poles, w;, ws are angular speeds 
respectively of the rotor and stator. 


4. CONTROL STRATEGY 
The active and reactive powers are written in the Park model as follows [15, 16]: 
Ps =V dsidstV qslqs 


: ‘ (8) 
Q s V qstds — V dslqs 


By the orientation of the stator flux along the axis d, we have: 
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@ > , Pas=%s (9) 
and 
Vas=Vs»>Vas=9 (10) 


The (8) becomes: 











Ps =V qsiqs (11) 
Q s —~V qslds 
By combining (5, 6, and 11), (11) becomes: 
P V 4 
= - 
Ss rig qr 
G27, 
s ae dr L. 
By substituting (6) for (5) and using (9), the rotor voltage can be rewritten as follows: 
; did . 
Var = RriartoL, a — 8050 Lyigr 
: 13 
| digr MV s = 
Var = Rrigr* OL + 89 Lyidr* 8 
t Ls 
M? Os 
Where: o=1-4"— ; g =l1-— 
LyLs Or 


5. FUZZY LOGIC CONTROL 
5.1. Principle of fuzzy logic method 

The principle of fuzzy control is based on the following three blocks: Fuzzification of input 
variables which are defined by linguistic variables (fuzzy sets or subsets), then the Inference where these 
variables are compared with predefined sets to determine appropriate response. And finally defuzzification, 
to convert the subset fuzzification in values, using the centroid defuzzification. The basic structure of our 
fuzzy controller is shown in Figure 2 [17, 18]. 






vit) 


Process 


Defuzzificatio 


Figure 2. Internal structure of a fuzzy controller 


With: R(t) is the reference signal, u(t) is the command signal, (t) is the output of the system to be controlled. 
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5.2. Regulating power by fuzzy logic 
The fuzzy controller is basically non-linear static input/output plot as show in Figure 3, 
the controller action can be written in the form: 


u=ke.e+kaeAe (14) 





Figure 3. Fuzzy logic regulator 


The fuzzy output is: 


y= ky U (15) 


Where: k, is the gain of the power error (Ps or Qs), ke is the gain of the variation of the power error (P; or Qs); 
e is the power error (P; or Q,), Ae is the variation of the error of (P; or Qs), u is the fuzzy output. 
The error, e is defined by: 


eae (16) 


e(k) = O5(k)-Os pep (K) 


Where: Ps ref is the reference stator power active, Qs ref is the reference reactive stator power, A@ is the 
variation of the error, it is approximated by: 


Ae(k) = e(k) — e(k -1) (17) 


5.3. Structure and application of a fuzzy adjustment of the stator power in a DFIG 
In this section, we supposed the fuzzy controller with five fuzzy subsets and constructed from the 
following condition rule: 


Rule (i): ifs is F',,thenuy is Fy, i=1,...,5 
Table 1 shows one of possible control rules based on five membership functions, where: NB is negative big, 
NS is negative small, ZE is zero, PS is positive small and PB is positive big. These acronyms are labels of 


fuzzy sets and their corresponding membership functions are depicted in Figure 4, [19]. To control the stator 
power in a DFIG we have using the block diagram of Figure 5 similar to that of Figure 3. 


Table 1.one of possible control rules based on five membership functions 








: Ae 
NB NS ZE PS PB 
NB ZE ZE PB PB PB 
NS ZE ZE PS PS PS 
ZE PS ZE ZE ZE NS 
PS NS NS NS ZE ZE 
PB NB NB NB ZE ZE 
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Figure 4. Membership functions of: (a) error e; (b) change of error Ae; (c) function of sortie u 





Generator PV 


Figure 5. Structure of a fuzzy adjustment of the stator power of a DFIG generator 


6. MODELING OF THE PHOTOVOLTAIC ENERGY CONVERSION SYSTEM 

This system consists of a PV solar generator directly connected with the chopper coverter through a 
bus continuous to feed the inverter converter with the DC voltage, which in turn feeds the rotor DFIG with the 
alternative voltage. A PV module consists of a number of solar cells connected in series and parallel to obtain 
the desired voltage and current output levels. As sunlight strikes a solar cell, the incident energy is converted 
directly into electrical energy. Single-diode mathematic model is applicable to simulate silicon photovoltaic 
cells, which consists of a photocurrent source, a nonlinear diode, and internal resistances, as shown 
in Figure 6. 





Figure 6. Equivalent diagram of a real photovoltaic cell [20] 
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IT pv=LTpna-lLa-TIp (18) 
e(V pvt pyRs) V pvt I pvRs 
T pv =T1 pn-~La =I pnr-To| exp 1 (19) 
py"? ye os mKT ¢ Rp 
—-E 
To =1.5*10>ex B (20) 
0 x10 (Z| 


Where: J, is the saturation current of the unlit junction; /, is the current of the cell; e is the charge of the 
electron (1.6 x10"°C); K is Boltzmann's constant (1.38 x10J / K); 7, is the junction temperature of the 
photovoltaic cell; m is the ideality factor of the junction; V,, is the voltage across the cell; R,; is the series 
resistance; R, is the shunt resistor. E,is the energy of the forbidden band for a semiconductor material gives, 
Eg is almost constant and Jp is therefore a function of the temperature. 


6.1. Modeling of the photovoltaic continuous bus 
The evolution of the photovoltaic DC bus voltage is given by the following equation: 


Los 
Vv pv — Jic,,a —) 
C pv 


The current in the photovoltaic DC bus comes from a node from which two currents flow: 


lo, = 1 pv ~thac (22) 


C pv 


Where: J.) is the current of the photovoltaic continuous bus; Jhac is the current of the chopper; Cy is the 
capacity of the photovoltaic continuous bus. 


6.2. Chopper booster 


The booster converter controls the electrical power in DC circuits with great flexibility and high 
efficiency. It consists of capacitors, inductance, and switches; all these devices consume no active power. 


: L Rina D 
“_,—SYYV_s\W\ > — 
iz 


i 
> 


x Vo 


i ple 


Figure 7. Electrical Boost converter Scheme [21] 





Vi 


During operation of the chopper, the transistor switches at a constant frequency fs with a closing time (aTs) 
and an opening time ((1-a) 7;.) where: 

- Ts is the switching period that equals (1/f) 

- a: the cyclic ratio of the switch (a € [0, 1]). 
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Figure 8. Periods of closing and opening a switch. 


Applying Kirchhoff's laws to the scheme of figure 7, we obtain the following systems of equations [22]: 
For the 1 period a. T; : 





pOa( i a) 2) 


dt 


jeg (0) = 62 POO ~ -56(0 (23) 


vp(t)= LEE = y,(1)-Rrin 


For the 2™ period (1-a).Ts: 


in (0) = = 1,(0)-12 0 
jest) =e. 2014 = 5, (0)-i0(0) 24) 
diz 


vi(t)=L~ == vi(t)— volt) Kriz 


To find a valid dynamic representation for the whole period 7;, we generally use the following expression: 


dx dx dx 
,.= Tt l-a )T (25) 
a Hert dilican ye 





Applying the relation (25) to the systems of equations (23) and (24), we find the Equations that govern 
the system over an entire period are obtained: 


dv; 
1 = T; = aT, (i; -i, )+(1-a@)T, (i; -iz ) 


d ee 
~ T, = -aT,ig + (1-@)T, (iz - ip) (26) 


Cc 





c2 





“Lr, = aT, (vj - Rziz )+(1-@)T, (vi — v9 - Rziz ) 


By arranging the terms of the preceding equations, (so that we can interconnect the Boost with the other 
simulation blocks), we obtain the dynamic modeling of the Boost converter: 





bo3 dvi(t) 
IL =1lj- Cc dt 
27) 
' dvolt ( 
io=(I-a)i, -02 v0) 
dit . 
Se +(l-a@ v9 + Rit 
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6.3. Modeling the chopper filter 
A filter is used to eliminate the harmonics generated by the chopper. The evolution of the filter 
voltage is given by the following equation [13]: 


di 
V fhac = R fhacl fhac + bie (28) 


Where: Rynac is the chopper filter resistance; Ljnac is chopper filter inductance. 


6.4. Chopper control 

Adjusting of the voltage according to the desired power (reference) is necessary to control the 
energy produced by the GPV photovoltaic generator with a modification of the variable control of the 
chopper (periodic ratio «) and servocontrol of the ifhac filter current (ifhac). In our study, the switch of the 
converter is controlled by a Pulse wave Modulation Large (PWM) signal, with a fixed frequency fs and a 
cyclic ratio a variable. The two operating phases of the switch S must be studied to extract the control law of 
the photovoltaic generator power output [23]: 





° Shac is closed (0<t<aT;): 
fae Tale” Yao” R phaci pee) (29 
° Siae is open (aT8<t<Ts): 
ai hae “Tpalrate* R maci rac) (30) 


Term (29) and (30) represent the equations of the chopper and its filter for a period a7s and (/-a)Ts 
respectively. So we have to find a valid approximated dynamic representation for both time slots. For this and 
from equation (25) we consider that the variation of the dynamic variable (imac) is of linear form, so the 
derivative of the dynamic variable (ifac) can be defined according to equation (25) over the two time periods 
aTs and (1-a)Ts: 


di a di a di hac 
l fhac ee l fhac aT + L fh (1-@ )r a 
a dt(ar,) dt((1-a)r ,) ; 


By replacing (29) and (30) by their value in relation (31) we obtain the equation that governs the system on the 
entire period: 


di fhac 1 


dt L fhac 





(av py-Vvde- R fraci prac) (32) 


The control law is given by the following expression and descripted in Figure 9: 


Vv 4] + Vdc 
pe fhac dc (33) 


V py 
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Eq. 


+ 
P oa PI —> MI oo» 
rep > va > a) ; 
- cc 


i hac j 
i lfhac Vac 





Figure 9. Chopper control. 


7. STRUCTURE OF THE MPPT COMMAND BY THE METHOD OF DISTURBANCE AND 
OBSERVATION « P & O» 

The "P & O" disturbance and observation method is a widely used approach in the search for 
Maximum Power Point Tracking (MPPT) because it is simple and requires only voltage and current 
measurements of the photovoltaic panels vpv and ipv respectively. It can be used to deduce the point of 
maximum power even during variations in illumination and temperature. 


7.1 Reference current calculation by program of the MPPT command 
Several algorithms have been proposed for the realization of the MPPT command for that we have 
chosen the following program: 


% Simple MPP "perturb and observe" tracking algorithm 
% using Boost DC-DC input current Iref as the control variable 
% Pold, Iref and Increment are initialized in Initialize MPPtrackIref.m 
% Input: power P to be maximized 

% Output: reference current 

function y = fen(Ppv) 

global Pold; 

global Iref; 

global Increment; 

IrefH = 5; % upper limit for the reference current 

IrefL = 0; % lower limit for the reference current 

Deltal = 0.02; % reference current increment 

P=Ppv 

if (P < Pold) 

Increment = -Increment; % change direction if P decreased 
end 

% increment current reference 
Iref=Iref+Increment*Deltal; 

% check for upper limit 

if (Iref > IrefH) 

Iref = IrefH; 

End 

% check for lower limit 

if (Iref < IrefL) 

Iref = IrefL; 

End 

% save power value 

Pold = P; 

% output current reference 

y = Iref; 


8. SIMULATIONS AND RESULTS 
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Simulation of the proposed system has been realized using Matlab/Simulink. The parameters of 
Simulation are given in appendix. The DFIG is connected directly through its stator to the network and his 
rotor to the cycloconverter (the inverter DC/AC) of the photovoltaic energy conversion system as show in 
Figure 10. The first part of simulation consists of the regulation system of the DC voltage obtained by 
the chopper. 









Electrical 
Nenworts 


Generator PV 


Figure 10. Control of the stator powers of the hybrid system by fuzzy logic 


In this section, the system described in Figure 10 was simulated in Matlab-Simulink. The first part of 
the simulation focuses on the regulation system of the DC voltage obtained by the chopper (regular at 330V) 
which feeds the inverter and then the rotor of the DFIG in along field changes of solar radiation. To achieve 
the control of the stator powers generated by the DFIG, the fuzzy logic controller is used as described in 
Figure 3. The second part focuses on the control system of the photovoltaic generator as described in Figure 6, 
where the goal is to extract the maximum power of solar sources. A part of this power is used to supply the 
rotor of the DFIG. The rest is injected into the grid. Finally we have a comparison between two systems, one is 
a classical system (rotor supply via electrical networks) and the other is a system that contains our proposal 
(rotor supply via a photovoltaic panel). 

The voltage of the chopper va is maintained at 330V as shown in Figure 11. A transition period be 
distinguished in the interval t = [0 0.04] s. After this period, the chopper voltage is maintained constant along 
the field changes of solar radiation as show in Figure 12. 








Insolation G (W/m2) 


DC Voltage vde (V) 





02 03 O4 OF 06 OF 08 09 
Time(s) 





Figure 11. The obtained chopper voltage vac Figure 12. Vareation of solar radiation G 
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In Figure 13, we have changes in the sequential reference current for solar radiation changes (Figure 
11). This current is used in the MPPT command program to track the maximum power point. The control of 
the active and reactive powers (P; , Qs) delivered by the generator DFIG is ensured by the fuzzy logic 
regulator. After simulating the system by means of Simulink / MATLAB, a fast acceptable variation of the 
voltage Vrape and the current [,as-¢ which supplies the rotor of the DFIG shown in Figures 14 and 15 has 
been found. 





reference current Iref (A) 














Time(s) 


Figure 13. Vareation of reference current of MPPT J,.r 
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Figure14. Rotor currents Figure 15. Rotor voltages 


Figure 16 and 17 show the regulation of active and active fixed power as a function of time at 
different times. At time t = 0.3s, the active and reactive power were adapted to the values of 4kW and Ikvar, 
respectively. A: t = [0.3 0.6] s: the active power was adapted to the value 5 kW, the reactive power to the value 
2kvar. At time t = 0.6s they decreased to the values of P;= 3 kW and Q, = 1.5kvar. 
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Figure 16. Stator-active power after control Figure 17. Stator-reactive power after control and 
and regulation regulation 
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For the energy produced by the photovoltaic generator show in Figure 18, it changes according to 
the vareation of solar radiation refer Figure 11 and reference current refer Figure 13. In parallel to 
the production of power output at the stator, we have consumption of internal energy by the rotor generator, 
this energy changes over time with change in energy levels required at the level of stator as show in 
Figure 19. At time t = 0.3s, the rotor power of DFIG consumed 317W. A: t = [0.3 0.6] s: this consumption up 
to 477W, and he down to 182W at time t = 0.6s. We can be used the residual energy produced by the 
photovoltaic generator after subtracting energy consumed from the rotor to charge the batteries and can be 
injected in to the grid as shown in Figure 20. 

In this section we have varation of the active and reactive powers stator Psa, Qsa and Ps, QO, of the 
two sytem, one is a classical system (rotor supply via electrical networks) and the other is a system that 
contains our proposal (rotor supply via a photovoltaic panel), respectively as show in Figure 21. A: t = [0 
0.25] s, we have tow systems work in the normal state. A: t = [0.25 0.4] s and t= [0.6 0.78] s, we have Lack 
in power of the electricity grid which causes a lack of production at the level of the classical system (P=4kw, 
QO.= \kvar and P,~=0kW, Qsa= Okvar). So the classical system is affected by network interruptions and our 
system proposed is remains stable. 











power photovotaic Ppv (W) 





Time(s) 


Figure 18. PV generated power 
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Figure 19. Rotor consumed power Figure 20. The rest of the power generated 
by the PV 
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Figure 21. Variation of active and reactive Stator power and reaction of tow systems deffirent after 
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9. CONCLUSION 
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In this paper, a hybrid energy conversion system combining photovoltaic and wind turbines as a 
scale alternative source of electrical energy is proposed. The power is generated by means of photovoltaic 
and wind energy while maintaining the energy efficiency and the reliability of the system. A complete 
description of the hybrid system was presented along with the detailed simulation results that prove its 
usefulness and In this context, maximum power was extracted under the operating conditions of the PV 
generator, and part of this energy was used to supply the DFIG rotor used in wind system and inject the rest 
of the energy into the grid. To adjust the power produced in the stator of DFIG was used the technique of the 
fuzzy logic because he is the simplest control technique Applied in practice. To extracted maximum power 
in photovoltaic systeme, we have used the tracking technique of maximum power point because she is simple 
and cost-effective. At the end, the entire circuit simulation was developed in a MATLAB/SIMULINK 
environment and the result based on the various simulations of the performance of system in different 
situations was presented and efficiency of the system proposed with classical system. 


APPENDIX 


Parameters of the system simulation 
Photovoltaic generator 

Insolation: G : [550 1000] W/m? 

Short circuit current: Js.: 5.45A 

Open circuit voltage: V,.: 22.2V 
Current at Pmax: J,.:4.95A 

Voltage at Pmax : V,: 17.2V 

Cell number in series : 2 

Cell number in Parallel : 44 


Wind generator 

Power: P: 7.5 KW 

Stator voltage: V;: 220V 

Rotor resistance: R;: 1.8Q 
Stator resistance: R,; : 0.95Q 
Rotor inductance: L,: 0.088H 
Stator inductance: L, : 0.094H 
Mutual inductance: M : 0.082H 
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